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The spectral function for an electron one-component plasma is calculated self-consistently using 
the GW' 9 ' approximation for the single-particle self-energy. In this way, correlation effects which 
go beyond the mean-field description of the plasma are contained, i.e. the collisional damping of 
single-particle states, the dynamical screening of the interaction and the appearance of collective 
plasma modes. Secondly, a novel non-perturbative analytic solution for the on-shell GW^ self- 
energy as a function of momentum is presented. It reproduces the numerical data for the spectral 
function with a relative error of less than 10% in the regime where the Debye screening parameter 
is smaller than the inverse Bohr radius, k < lttj 1 . In the limit of low density, the non-perturbative 
self-energy behaves as n 1 ^ 4 , whereas a perturbation expansion leads to the unphysical result of a 
density independent self-energy [W. Fennel and H. P. Wilfer, Ann. Phys. Lpz. 32, 265 (1974)]. 
The derived expression will greatly facilitate the calculation of observables in correlated plasmas 
(transport properties, equation of state) that need the spectral function as an input quantity. This 
is demonstrated for the shift of the chemical potential, which is computed from the analytical 
^ ■ formulae and compared to the GW '"■'-result. At a plasma temperature of 100 eV and densities 

O ' below 10 21 cm" 3 , both approaches deviate less than 10% from each other. 
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INTRODUCTION 



The many-particle Green function approach 1] allows for a systematic study of macroscopic properties of correlated 
systems. Green functions know a long history of applications in solid state theory [2j, nuclear [3J, and hadron physics 
O ,| and also in the theory of strongly coupled plasmas [f|. In the latter case, optical and dielectric properties [0,01 have 

^ . been studied using the Green function approach, as well as transport properties like conductivity Q and stopping 
power 0, [l(| , and the equation of state [ll[ . Modifications of these quantities due to the interaction among the 
constituents can be accessed, starting from a common starting point, namely the Hamiltonian of the system. 
^-j 1 The key quantity to electronic properties in a correlated many-body environment is the electron spectral function 
' ^4.(Pj i- c - the probability density to find an electron at energy (frequency) ui for a given momentum p. It is related 
to the retarded electron self-energy E(p, u> + iS), S > via Dyson's equation 

rsl ! nr \ -2Im£(p,w + i5) 

[u - e P - Re £(p, lu)} 2 + [Im E(p, w + i5)f 

■ Here, the single-particle energy 

<~1: e P =P 2 -/ic (2) 

OO ' 

has been introduced, /j, e is the electron chemical potential. Note that here and throughout the paper, the Rydberg 
system of units is used, where % = 1, m c = 1/2, and e 2 /47reo = 2. Furthermore, the Boltzmann constant fce is set 
' equal to 1, i.e. temperatures are measured in units of energy. 
J> , The self-energy describes the influence of correlations on the behaviour of the electrons. A finite, frequency depen- 
dent self-energy leads to a finite life-time of single-particle states and a modification of the single-particle dispersion 
relation. Hence, the calculation of the electron self-energy is the central task if one wants to determine electronic 

■ properties, e.g. those mentioned above. 

The Hartree-Fock approximation lT2j represents the lowest order in a perturbative expansion of the self-energy in 
terms of the interaction potential [131 ]. Being a mean- field approximation, effects due to correlations in the system 
cannot be described. Examples are the appearance of collective modes, the energy transfer during particle collisions, 
and the quasi-particle damping. The next order term is the Born approximation, where binary collision are taken 
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into account via a bare Coulomb potential. However, the Born approximation leads to a divergent integral, due to 
the long-range Coulomb interaction. Therefore, the perturbation expansion of the self-energy has to be replaced by 
a non-perturbative approach, accounting for the dynamical screening of the interparticle potential. 

A non-perturbative approach to the many-particle problem is given by the theory of Dyson [14j and Schwinger 
[IM li^t generalized to finite temperature and finite density [l7| . An excellent introduction to Dyson-Schwinger 
equations can also be found in j4j. The Dyson-Schwinger equation for the self-energy £ contains the full Green 
function G, the screened interaction W and the proper vertex T. Since each of these functions obeys a different 
Dyson-Schwinger equation itself, involving higher order correlation functions, the Dyson-Schwinger approach leads to 
a hierarchy of coupled integro-differential equations. In order to provide soluble equations, this hierarchy has to be 
closed at some level, i.e. correlation functions of a certain order either have to be parametrized or neglected. 

One such closure of the Dyson-Schwinger hierarchy consists in neglecting the vertex, i.e. the three-point function, 
and considering only the particle propagators and their respective self-energies, i.e. two-point functions. One arrives at 
the so-called GVF-approximation, introduced in solid-state physics by Hedin [1S| . Hedin was led by the idea, to include 
correlations in the self-energy by replacing the Coulomb potential in the Hartree-Fock self-energy by the dynamically 
screened interaction W. In this way, one obtains a self-consistent, closed set of equations for the self-energy, the 
polarization function II, the Green function and the screened interaction. 

It can be shown [l9j] , that the GW approximation belongs to the so-called ^-derivable approximations [13, HH • As 
such, it leads to energy, momentum, and particle number conserving expressions for higher order correlation functions. 
It was successfully applied in virtually all branches of solid state physics. An overview on theoretical foundations and 
applications of the GW approximation can be found in the review articles [HI, [H, [Hj] • 

The drawback of the GW approximation is, that the Ward-Takahashi identities are violated. Ward-Takahashi 
identities provide an exact relation between the vertex function T, i.e. the effective electron-photon coupling in the 
medium, and the self-energy and follow from the Dyson-Schwinger equations. They reflect the gauge invariance of the 
theory. In GW, they are violated simply because corrections to the vertex beyond zero order are neglected altogether. 
This issue touches on a fundamental problem in many-body theory and field theory, namely the question how to 
preserve gauge invariance in an effective, i.e. approximative theory, without violating basic conservation laws. A 
detailed analysis of this question with application to nuclear physics is presented in a series of papers by van Hees 
and Knoll [ll, [11, S3 • 

Approximations for the self-energy, that also contain the vertex are often referred to as GWT approximations. 
An example can be found in Ref. [28j . where the spectral function of electrons in aluminum is calculated using a 
parametrized vertex function. An interesting result obtained in that work is that vertex corrections and self-energy 
corrections entering the polarization function, largely cancel. This can be understood as a consequence of Ward- 
Takahashi identities. Thus, and in order to reduce the numerical cost, it is a sensible choice to neglect vertex 
corrections altogether, and to keep the polarization function on the lowest level, i.e. the random phase approximation 
(RPA) which is the convolution product of two non-interacting Green functions in frequency-momentum space. The 
corresponding self-energy is named the GW^ self-energy and has been introduced by von Barth and Holm [29| . 
who were also the first to study the fully self-consistent GW approximation J3C|. Throughout this work, the GW^ 
self-energy will be analyzed. 

Having been used in solid state physics traditionally, the GW /< -°- ) -method was also applied to study correlations in 
hot and dense plasmas, recently. The equation of state [UUll, as wei l as optical properties of electron hole plasmas 
in highly excited semiconductors [HI], and dense hydrogen plasmas 0] were investigated. 

In general, the calculation of such macroscopic observables of many-particle system involves numerical operations 
that need the spectral function as an input. Since the self-consistent calculation of the self-energy, even in GW^- 
approximation, is itself already a numerically demanding task, it is worth looking for an analytic solution of the GW^ 
equations, which reproduces the numerical solution at least in a certain range of plasma parameters. Such an analytic 
expression then also allows to study the self-energy in various limiting cases, such as the low density limit or the limit 
of high momenta, which are difficult to access in the numerical treatment. Furthermore, an analytic expression, being 
already close to the numerical solution permits the calculation of the full GW^ self-energy using only few iterations. 

Anal ytic al expressions for the single-particle self-energy have already been given by other authors, e.g. Fennel and 
Wilfer |34j and Kraeft et al. [HI]. They calculated the self-energy in first order of the perturbation expansion with 
respect to the dynamically screened potential. Besides being far from the converged GW^ self-energy, their result 
is independent of density, i.e. the single-particle life-time is finite even in vacuum. As shown in [35j |. this unphysical 
behaviour is a direct consequence of the perturbative treatment. By using a non-perturbative ansatz, an expression 
for the self-consistent self-energy in a classical one-component plasma was presented, that reproduces the full GW^ 
self-energy at small momenta, i.e. for slow particles. The behaviour of the quasi-particle damping at larger momenta 
remained open and will be investigated in the present work. Secondly, based on the information gathered about the 
low and high momentum behaviour, an interpolation formula will be derived, that gives the quasi-particle damping 
at arbitrary momenta. 



The work is organized in the following way: After a brief outline of the GW^-approximation in the next section, 
numerical results will be given in section IIIII for the single-particle spectral function for various sets of parameters 
electron density n e and electron temperature T. In section ITVl the analytic expression for the quasi-particle damping 
width is presented and comparison to the numerical results are given. Section [V] deals with the application of the 
derived formulae to the calculation of the chemical potential as a function of density and temperature. An appendix 
contains the detailed derivation of the analytic self-energy. As a model system, we focus on the electron one-component 
plasma, ions are treated as a homogeneously distributed background of positive charges (jellium model). 



II. SPECTRAL FUNCTION AND SELF-ENERGY 



We start our discussion with the integral equation for the imaginary part of the single-particle self-energy in GW^ 
approximation: 

Im£(p,w + i<S) = — j^rY" / —V(q)A(p-q,uj-uj')lme^ >A ((i,uj')n B (uj')n F (uj-uj'). (3) 
™fM q J-oc 27r 

V(q) = 8ir/q 2 flo is the Fourier transform of the Coulomb potential with the normalization volume f2o- It is multiplied 
by the inverse dielectric function in RPA, 

I \ 1 \ \ " "F( £ P-q/2) ~ ^F(£p+q/2) , . 

e R ,PA(g,w) = 1 - V(q) > ■ . (4) 

- U + £p-q/2 - £p+q/2 

to account for dynamical screening of the interaction. Furthermore, the Fermi-Dirac and the Bose-Einstein distribution 
function, rip/B^) — [exp(e//cBT) ±1] 1 , were introduced. Note, that the dielectric function is only determined once, 
at the beginning of the calculation. In particular, the single-particle energies £ p = p 2 — /i c entering equation Q are 
determined from the non-interacting chemical potential, whereas during the course of the self-consistent calculation 
of the self-energy, the chemical potential is recalculated at each step via inversion of the density relation 



(5) 



Using the self-consistent chemical potential also in the RPA polarization function leads to violation of the /-sum rule, 
i.e. conservation of the number of particles. 

The real part of the self-energy is obtained via the Kramers-Kronig relation [2j All quantities (spectral function, 
self-energy, and chemical potential) have to be determined in a self-consistent way. This is usually achieved by solving 
equations (jTJ) , , and iteratively.The numerical algorithm is discussed in detail in Ref. [351 ]. 



III. NUMERICAL RESULTS 



The spectral function was calculated for the case of a hot one-component electron plasma. Temperature and density 
were chosen such, that the plasma degeneracy parameter 

*=£• < 6 » 

is always larger than 1, i.e. the plasma is non-degenerate. Furthermore, the temperature is fixed above the ionization 
energy of hydrogen, T » 1 Ry, such that bound states can be neglected in the calculations. At lower temperature, 
bound states have to be include, e.g. via the t-matrix. For an application in electron-hole plasmas, see reference [33l ]. 
The electron density is adjusted such that the the plasma coupling parameter 

2 /47m e x 1/3 



which gives the mean Coulomb interaction energy compared to the thermal energy, is smaller than 1 in all calculations, 
i.e. we are in the limit of weak coupling. 



In figure [TJ we show contour plots of the spectral function in frequency and momentum space for two different 
densities (n c — 7 x 10 21 cm -3 , upper graph) and (n e — 7 x 10 25 cm~ 3 , lower graph). The temperature is set to 
T = 1000 eV in both calculations. The free particle dispersion oj = e p is shown as solid black line. 

In the first case, the plasma is classical (0 — 7.5 x 10 2 ) and weakly coupled (r = 4.4 x 10~ 3 ). The spectral function is 
symmetrically broadened and the maximum is found at the free dispersion, i.e. there is no notable quasi-particlc shift 
at the present conditions. At increasing momentum, the width of the spectral function decreases, and the maximum 
value increases; the norm is preserved. 

The situation changes, when we go to higher densities, cf. the lower graph in figure [T] The chosen parameters are 
typical solar core parameters (36|. The degeneracy parameter is now — 1.6 and the coupling parameter is T — 0.096. 
The increased degeneracy and coupling result in a significant modification of the spectral function as compared to 
the low density case: A shift of the spectral function's maximum to smaller frequencies is observed, the Hartree-Fock 
shift. The shift due to dynamic correlations is still small at the present conditions; it becomes important in strongly 
degenerate systems [29) . Furthermore, shoulders appearing in the wings of the main quasi-particle peak at small 
momenta, indicate the excitation of new quasi-particles, so-called plasmarons [37j. They can be seen more clearly in 
figure [7] (solid curve). The plasmaron satellites are separated from the main peak by roughly the plasma frequency 
utpi = 4^/7m c which is about 23 Ry at the present density. In the former case of lower density no plasmarons appear, 
only a featureless, single resonance is obtained. At higher momenta, the plasmarons merge into the central peak. As 
in the low-density case, the position of the maximum approaches the single-particle dispersion, due to the decreasing 
Hartree-Fock shift at high p. Again, the width of the spectral function decreases with increasing momentum. 
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FIG. 1: (Color online) Contour plots of the spectral function as a function of momentum and frequency. The colour scale is 
logarithmic. Results are shown for two different densities (n — 7 x 10 21 cm -3 , upper graph) and (n = 7 x 10 25 cm -3 , lower 
graph). For these parameters, the plasma coupling parameter is T — 4.4 x 10 -3 , and F = 9.6 x 10~ 2 , respectively. The 
degeneracy parameter is 8 — 7.5 x 10 2 and 6 = 1.6, respectively. The black line indicates the free particle dispersion uj = s p . 

This is visible more clearly in figures [2] - [71 Here, the solid curves represent the GW^ spectral function as a 
function of frequency. Results are shown for three different momenta, p = Oa^ 1 (a), p — 50 a^ 1 (b), and p = 100 a^ 1 



(c). Two different temperatures are considered, T = fOOeV (figures and T = fOOOeV (figures [M?]) and for each 
temperature three different densities are studied. With increasing momentum p, the spectral function becomes more 
and more narrow, converging eventually into a narrow on-shell resonance, located at the unperturbed single-particle 
dispersion uj + /j, c = p 2 . 

As a general feature, one can observe an increase of the spectral function's width with increasing density and 
with increasing temperature. The increase with density is due to the increased coupling, while the increase with 
temperature reflects the thermal broadening of the momentum distribution function n-p{uj) that enters the self-energy 
and thereby also the spectral function. From these results, we see that the spectral function has a quite simple form 
in the limit of low coupling, i.e. at low densities and high temperatures. 

The numerical results are compared to a Gaussian ansatz for the spectral function, shown as the dashed curve in 
figures [2] - [3 The explicit form of the Gaussian is given as equation (fTTj) . below. It's sole free parameter is the width, 
denoted by a p . An analytic expression for a p will be derived in section HVl The coincidence is in general good at high 
momenta, whereas at low momenta, the spectral function deviates from the Gaussian. In particular, the steep wings 
and the smoother plateau that form at low momenta are not reproduced by the Gaussian. Also, the plasmaron peaks 
appearing in the spectral function at high density (see figure cannot be described by the single Gaussian. 

Determining a p via least-squares fitting of the Gaussian ansatz to the numerical data at each p leads to the solid 
curve in figure [5]), obtained in the case of n = 7 x 10 20 cm~ 3 and T = 100 eV. Starting at some finite value Co at 
p = 0, the width falls off slowly towards larger p. The dashed curve shows a p as obtained from the analytic formula 
that will be derived in the now following section. 



IV. ANALYTICAL EXPRESSION FOR THE QUASI-PARTICLE SELF-ENERGY 

The solution of the GW^- equation ^ requires a considerable numerical effort. So far (see e.g. the work by 
Fennel and Wilfer in [Uj]), attempts to solve the integral ([3|) analytically were led by the idea to replace the spectral 
function on the r.h.s. by its non-interacting counterpart, A^ (p, uj) — 2TrS(e p — uj), i.e. to go back to the perturbation 
expansion of the self-energy and neglect the implied self-consistency. At the same time, the inverse dielectric function 
is usually replaced by a simplified expression, e.g. the Born approximation or the plasmon-pole approximation [lzl ]. 
Whereas the second simplification is indispensable due to the complicated structure of the inverse dielectric function, 
the first one, i.e. the recursion to the quasi-particle picture, is not necessary, as was shown by the author in Ref. [35|. 
In fact, the result that one obtains in the quasi-particle approximation is far from the converged result, at least in 
the high temperature case. Secondly, using the quasi-particle approximation, the imaginary part of the self-energy is 
not density dependent, i.e. a finite life-time of the particle states is obtained even in vacuum. This unphysical result 
can only be overcome if one sticks to the self-consistency of the self-energy, i.e. if one leaves the imaginary part of 
the self-energy entering the r.h.s. of equation (J3J) finite. 

Using the statically screened Born approximation, which describes the binary collision among electrons via a 
statically screened potential, a scaling law Im S(p, o>Q p (p)) oc T 3 / 4 was found [35]. Hence, the spectral functions 
width vanishes when the plasma coupling parameter T (see equation (JT])) tends to 0. An expression for the self-energy 
was found, that reproduces the converged GW^ calculations at small momenta, p <C K. At higher momenta, the 
derived expression ceases to be valid. 

In this work, a different approximation to the dielectric function is studied, namely the plasmon-pole approximation 
This means, that the inverse dielectric function is replaced by a sum of two delta-functions that describe the 
location of the plasmon resonances, 

Ime RPA(q) w ') — > Imep^ A (q,u/) =---£- [S(uj - ujq) + 5(uj + Uq)] . (8) 

Z CiJ q 

For classical plasmas, the plasmon dispersion u q can be approximated by the Bohm-Gross dispersion relation (38j [. 

^^(l + jj^+tf*. (9) 

Many-particle and quantum effects on the plasmon dispersion have recently been studied in [39j . 

The plasmon-pole approximation allows to perform the frequency integration in equation ([3]), resulting in the 
expression 



UJ a 



A(p - q, uj - ujq)n B (w q ) exp(w q /T) - A(p - q, uj + w q )n B (~w q ) exp(-cj q /T) 



(10) 



We will first study the case of high momenta, i.e. momenta that are large against any other momentum scale or 
inverse length scale, such as the mean momentum with respect to the Boltzmann distribution, p = y/3T/2 or the 
inverse screening length k = yJ&irn/T. 

As discussed in the previous section, the numerical results for the spectral function at high momenta can well be 
reproduced by a Gaussian. Thus, we make the following ansatz for the spectral function: 

A Gauss (p,^) = -^— exp(-^ . (11) 

Note, that only the Hartree-Fock contribution to the real part of the self-energy appears, the frequency dependent 
part is usually small near the quasi-particle dispersion, 

W Q p (p) = e p + Re E(p, w) L =wQ p (p) , (12) 

which therefore can be approximated as ui Qp (p) = e p + Re£ HF (p). In the following, we make use of the knowledge 
about the width parameter a p that we gathered already through simple least-squares fitting of the Gaussian ansatz 
to the spectral function in order to solve the integrals in equation (fTU|) . 

First, we replace the spectral function on the r.h.s by the Gaussian ansatz (fTTj) and evaluate the emerging equation 
at the single-particle dispersion w^ p (p). By claiming that the Gaussian and the spectral function have the same 
value at the quasi-particle energy, we identify a p — ^/7r/2ImI](p,ti; QP (p))). Figure [H] shows, that the quasi-particle 
damping a p is a smooth function of p, that varies only little on the scale of the screening parameter k. Since the 
latter defines the scale on which contributions to the g-integral are most important, we can neglect the momentum 
shift in the self-energy on the r.h.s., i.e. we can replace the spectral function on the r.h.s of equation (fTUf as 

A(p - q, e p + E HF (p) ± o, q ) -> exp (- (£p ± = £p - q)2 ) , (13) 

<j p V 2a p J 

and can now perform the integral over the angle 9 between the momenta p and q, 



d cos 9 exp 



(e p ± w q - p 2 - q 2 + 2pqcos9 + /x e ) 5 



2a 2 



TT 

2p~q 



Erf 



(p + q f -p 2 T^ q 
V2a p 



Erf 



(p - q) 2 ~p 2 ±Uq 
V2a p 



(14) 



The remaining integration over the modulus of the transfer momentum q can be performed after some further 
approximations, explained in detail in appendix El For large p, one finally obtains the transcendent equation 



I-'^m/^ [n B (w p i) exp(w p i/T) -n B (-w p i) exp(-w pl /T)] 



7T T 

22~p 



Mk 2 p 2 I<j 2 p ) 



The solution of this equation can be expanded for large arguments of the logarithm, yielding 



ttT 
2p 
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<f(p) > <p(p) 



£(p)-ln£(p) 



lng(p) lng(p) ln|(p) 31n 2 C(p) ln 2 ^(p) In 3 g P ) 



£(p) = ln(W -K P 2 exp{A/T)/T) , A = -1.3357^ [n B (wpi) exp(o; pl /T) - n B {-u pl ) exp(-Wpi/T)] 

7T Z 



(15) 



o( P y 



(16) 



Equation (fT6|) is a solution of equation p5|) provided the argument of the inner logarithm is larger than Euler's 
constant e, i.e. \J~^Kp 2 exp(A/T)/T > e, i.e. at large p. The case of small p, where the previous inequality does not 

hold, has to be treated separately, see appendix [Bl 

Together with an expression for the quasi-particle damping at vanishing momentum taken from [35| and scaled such 
that the maximum of the spectral function at p = is reproduced, ctq = — 7rV kT/2, an interpolation formula (Pade 



formula) was derived that covers the complete p-range: 



Padc 



£(p) - ln(e + 



Qq + aiP 
1 + bip + hip 1 

2 



<f(p) > <P(P) = 



£(p)-lnf(p) 



ln£(p) ln|(p) , ln£(p) 31n 2 £(p) , In 2 |(p) , In 3 |(p) 



£(p) ^ 2 (P) £(p) 2£3( p ) 2C 2 (p) 3f3( p ) 



-Kp 2 exp(yl/T)/r) 



ai = — k 



(I) 



3/2 
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- 2T 



(17) 



The function £(p) in the last equation differs from £(p) in equation (|16p in that Euler's constant e ~ 2.7183 has been 
added to the argument of the logarithm. In this way, the function (p{p) is regularized at small p and tends to 1 at 
p = 0, i.e. the quasi-particle damping goes to the correct low-p limit. At large p, this modification is insignificant, 
since the original argument rises as p 2 . For the detailed derivation, see appendix [Bl 

Expression (|17[l . used in the Gaussian ansatz (fTT|) . leads to a spectral function that well reproduces the numerical 
data from full GVF^-calculations: Figure [8] (dashed curve) shows the effective quasi-particle damping width a p as a 
function of momentum p for the case n — 7 x 10 20 cm -3 and T — 100 eV. The solid curve gives the best-fit value for 
cr p obtained via least-square fitting of the full GW^-calculations assuming the Gaussian form (111|) . see section HTT1 
Both curves coincide to a large extent. The largest deviations are observed in the range of p ~ 20 a^ 1 . At this point, 
the validity of expression (fTfj|) as the solution of equation (fT5"|) ceases, since the argument of the logarithm becomes 
smaller than e. As already mentioned, we circumvented this problem by regularizing the logarithms, adding e to its 
argument. The deviation at p ~ 20 a^ 1 of up to 15% is a residue of this procedure. At higher momenta, the deviation 
is generally smaller than 10% and the analytic formula evolves parallel to the fit parameters. 
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FIG. 2: Spectral function in GW^-approximation (solid lines) and Gaussian ansatz (dashed lines) with quasi-particle damping 
width Up taken from equation (|17|l for three different momenta. Plasma parameters: n = 7 x 10 19 cm~ 3 , T — 100 eV. The 
plasma coupling parameter is T — 1.0 x 10~ 2 , the degeneracy parameter is 6 = 1.6 x 10 3 , the Debye screening parameter is 
K = 6.0 x 10~ 3 a" 1 . 



At smaller densities, the correspondence is even better as can be seen by comparing the spectral functions shown 
in figures [5] - [7] The dashed curves give the Gaussian ansatz for the spectral function with the quasi-particle width 
taken from the interpolation formula (fT7|) . As a general result, the analytic expression for the quasi-particle damping 
Up leads to a spectral function that nicely fits the numerical solution for the spectral function at least at finite p. At 
very small values of p, the overall correspondence is still fair, i.e. the position of the maximum and the overall width 
match, but the detailed behaviour does not coincide. In particular, the steep wings and the central plateau, that 
forms in the GIF' ' -calculation, is not reproduced by the one-parameter Gaussian. For this situation, the analytic 
formula for self-energy given in (3f| should be used instead. 

By comparing the numerical data for the spectral function to the Gaussian ansatz at different densities, it is found, 
that the Gaussian spectral function is a good approximation as long as the Debye screening parameter k is smaller 
than the inverse Bohr radius, K < lajj 1 - This becomes immanent by comparing figures [H and [7J In the first case 
(ri c = 7 x 10 23 cm- 3 ,T = 1000 eV), we have K = 0.19, while in the second case (n e = 7 x 10 25 cm- 3 ,T = 1000 eV), 
K = 1.9 is found. As already noted in the discussion of the numerical results in section [TTl] in the case of increased 
density, the plasmaron satellites appear as separate structures in the wings of the central quasi-particle peak, whereas 
they are hidden in the central peak at smaller densities. Therefore, a single Gaussian is not sufficient to fit the 
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FIG. 3: Spectral function in GW^-approximation (solid lines) and Gaussian ansatz (dashed lines) with quasi-particle damping 
width a p taken from equation (|17|1 for three different momenta. Plasma parameters: n = 7 x 10 20 cm -3 , T = 100 eV. The 
plasma coupling parameter is T — 2.1 x 10~ 2 , the degeneracy parameter is 6 = 3.5 x 10 2 , the Debye screening parameter is 
k — 1.9 x lQ-'d^ 1 . 
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FIG. 4: Spectral function in GW^-approximation (solid lines) and Gaussian ansatz (dashed lines) with quasi-particle damping 
width Up taken from equation (|17|l for three different momenta. Plasma parameters: n = 7 x 10 21 cm -3 , T — 100 eV. The 
plasma coupling parameter is T — 4.4 x 10~ 2 , the degeneracy parameter is 6 = 7.5 x 10 1 , the Debye screening parameter is 
K = 6.0 x 10 _2 a" 1 . 



spectral function at increased densities. Since the position of the plasmaron peak is given approximatively by the 
plasma frequency uj p \, whereas the width of the central peak at small p is just the quasi-particle width ctoj we can 
identify the ratio of these two quantities, — Lo p \/ao cx y/~K as the parameter which tells us if plasmaron peaks appear 

1 /2 

separately [u p i > —<Jo) or not (w p i < — erg). Since the plasma frequency increases as a function of n e , whereas the 
quasi-particle width scales as n 1 / 4 (c.f. equation (|17p ). the transition from the single peak behaviour to the more 
complex behaviour including plasmaron resonances, appears at increased density. Neglecting numerical constants of 
order 1 in the ratio of plasma frequency to damping width, we see that — Wpi/cn < 1 is equivalent to k < 1, which 
was our observation from the numerical results. Therefore, we can identify the range of validity of the presented 
expressions for the spectral function and the quasi-particle damping. It is valid for those plasmas, where we have 
densities and temperatures, such that k < 1. 

The physical origin of the requirement n < 1 can be understood in the following way [35| . At length scales smaller 
than the Bohr radius, one typically expects quantum effects, e.g. Pauli blocking. These effects are not accounted for 
in the derivation of the quasi-particle damping. Therefore, it appears to be a logical consequence that the validity of 
the results is limited by the length scale at which typical quantum phenomena become important. 

The regime of validity of the analytic formula can also be expressed via the plasma coupling parameter and the 
temperature as T < T~ 2 / 3 . Since we restrict ourselves to plasma temperatures, where bound states can be excluded, 
i.e. T > 1 Ry, this is equivalent to saying T < 1. 

Although the correspondence between the accurate GW^ calculations and the parametrized spectral function at 
small momenta is not as good as in the case of large momenta, the parametrized spectral function can be applied in 
the regime of validity to the calculation of plasma observables without introducing too large errors. As an example, 
this will be shown for the case of the chemical potential /i in the next section. 
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FIG. 5: Spectral function in GW^'-approximation (solid lines) and Gaussian ansatz (dashed lines) with quasi-particle damping 
width a p taken from equation (|17p for three different momenta. Plasma parameters: n = 7 x 10 21 cm~ 3 , T = 1000 eV. The 
plasma coupling parameter is T — 4.4 x 10~ 3 , the degeneracy parameter is 8 — 7.5 x 10 2 , the Debye screening parameter is 
K = 1.9 X 10" 2 a" 1 . 
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FIG. 6: Spectral function in GW^'-approximation (solid lines) and Gaussian ansatz (dashed lines) with quasi-particle damping 
width (j p taken from equation (|17p for three different momenta. Plasma parameters: n = 7 x 10 23 cm -3 , T = 1000 eV. The 
plasma coupling parameter is T — 2.1 x 10 -2 , the degeneracy parameter is 8 = 3.5 x 10 1 , the Debye screening parameter is 

K ~ 1.9 X 10 _1 flg 1 . 



V. APPLICATION: SHIFT OF THE CHEMICAL POTENTIAL 



To demonstrate the applicability of the presented formulae for quick and reliable calculations of plasma properties, 
we calculate the shift of the electron's chemical potential A/i = \i — /if reo , i.e. the deviation of the chemical potential of 
the interacting plasma [i from the value of the non-interacting system ^f rcc . The chemical potential of the interacting 
system [i is obtained by inversion of the density as a function of T and fj,, equation ([5|). The free chemical potential 
/ifree is obtained in a similar way by inversion of the free density 

n hcc (T, /Xfree) = 2^n F (e p ~ /Ltf re e) . (18) 
P 

Figure[9]shows the shift of the chemical potential as a function of the plasma density n for a fixed plasma temperature 
T = 100 eV. Results obtained by inversion of equation (JSJ using the parametrized spectral function (111|) with the 
quasi-particle damping width taken from equation (| 1 T[) (solid curve) are compared to those results taking the numerical 
GW(°) spectral function (dashed curve). 

The GW<°) result gives slightly smaller shifts than the parametrized spectral function, i.e. the usage of the analytical 
damping width leads to an overestimation of the the shift of the chemical potential. However, the deviation remains 
smaller than 20% over the range of densities considered here, i.e. for n < 1. At small densities, i.e. for n < 10 20 cm~ 3 , 
the parametrized spectral function yields the same result as the full GW^ calculation. 

The deviation at increased density can be reduced by improving the parametrization of the spectral function at 
small momenta. To this end, the behaviour of the quasi-particle damping width at high momenta, equation (I16|) 
should be combined with the frequency dependent solution for a p at vanishing momentum, as presented in Ref. [35| |. 
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FIG. 7: Spectral function in GW^° '-approximation (solid lines) and Gaussian ansatz (dashed lines) with quasi-particle damping 
width a p taken from equation (|17p for three different momenta. Plasma parameters: n = 7 x 10 25 cm -3 , T = 1000 eV. The 
plasma coupling parameter is F = 9.6 x 10 -2 , the degeneracy parameter is 8 — 1.6, the Debye screening parameter is n — 1.9 <ig 1 . 
Here, the Gaussian fit is no longer sufficient due to the appearance of plasmaron resonances in the spectral function (shoulders 
at lo ~ -30 Ry and u ~ 20 Ry). 
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FIG. 8: Effective quasi-particle damping a p as a function of momentum p for plasma density n — 7 x 10 20 cm -3 and temperature 
T = 100 eV. The fit-parameters for the Gaussian fit to the full GW^-calculations are given as dots, the solid line denotes the 
analytic interpolation formula (|17[1 . 

However, this task goes beyond the scope of this paper, where we wish to present comparatively simple analytic 
expressions for the damping width that yield the correct low density behaviour of plasma properties. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper, the GW-^ ' -approximation for the single-particle self-energy was evaluated for the case of a classical 
one-component electron plasma, with ions treated as a homogeneous charge background. A systematic behaviour of 
the spectral function was found, i.e. a symmetrically broadened structure at low momenta and convergence to a sharp 
quasi-particle resonance at high p. At increased densities, plasmaron satellites show up in the spectral function as 
satellites besides the main peak. 

In the second part, an analytic formula for the imaginary part of the self-energy at the quasi-particle dispersion 
w^ p (p) = £ p +E HF (p) was derived as a two-point Pade formula that interpolates between the exactly known behaviour 
at p = and p — > oo. The former case was studied in (35|, while an expression for the asymptotic case p — > oo was 
derived here. The result is summarized in equation (|17|) . In contrast to previously known expressions for the quasi- 
particle damping, based on a perturbative approach to the self-energy (34j , the result presented here shows a physically 
intuitive behaviour in the limit of low densities, i.e. it vanishes when the system becomes dilute. Using the Gaussian 
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FIG. 9: Shift of the chemical potential as a function of the plasma density for a plasma temperature T = 100 eV. Results 
for Afi using the parametrized spectral function (solid line) are compared to full numerical calculations, using the GW^- 
approximation. 

ansatz (fTT|) for the spectral function in combination with the quasi-particle width leads to a very good agreement 
with the numerical data for the spectral function in the range of plasma parameters, where k < lag 1 ; the relative 
deviation is smaller than 10% under this constraint. 

Thus, a simple expression for the damping width of electrons in a classical plasma has been found, that can be 
used to approximate the full spectral function to high accuracy. This achievement greatly facilitates the calculation 
of observables that take the spectral function or the self-energy as an input, such as optical properties (inverse 
bremsstrahlung absorption), conductivity, or the stopping power. 

Furthermore, it was demonstrated that the derived expressions allow for quick and reliable calculations of plasma 
properties without having to resort to the full self-consistent solution of the GW^ '-approximation. As an example, 
the shift of the chemical potential was calculated using the parametrized spectral function and compared to GW^ 
results. For densities of n < 10 21 cm~ 3 , both approaches coincide with a relative deviation of less than 10%, going 
eventually up to 20% as the density approaches 10 22 cm~ 3 . At low densities both approaches give identical results. 
This shows the extreme usefulness of the presented approach for the calculation of observables via the parametrized 
spectral function. 

As a furhter important application of the results presented in this paper, we would like to mention the calculation of 
radiative energy loss of particles traversing a dense medium, i.e. bremsstrahlung. A many-body theoretical approach 
to this scenario is given by Knoll and Voskresensky (40l | , using non-equilibrium Green functions. They showed, that a 
finite spectral width of the emitting particles leads to a decrease in the bremsstrahlung emission. This effect is known 
as the Landau-Pomeranchuk-Migdal effect [4l|, [42j. It has been experimentally confirmed in relativistic electron 
scattering experiments using dense targets, e.g. lead [H, HI]. I n Hi it is shown that also thermal bremsstrahlung 
from a plasma is reduced due to the finite spectral width of the electrons in the plasma. In the cited papers, the quasi- 
particle damping width was either set as a momentum- and energy independant parameter (in |40(), or calculated 
self-consistently using simplified approximations of the GW^ theory (in [45J), which itself is a very time-consuming 
task and prohibited investigations over a broad range of plasmas parameters. Now, based on this work's results, 
calculations on the level of full GW^ become feasible, since analytic formulae have been found that reproduce the 
GW^ self-energy. Effects of dynamical correlations on the bremsstrahlung spectrum can be studied starting from a 
consistent single-particle description via the GW^ ' self-energy. 
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APPENDIX A: ANALYTIC SOLUTION FOR THE GW (()} SELF-ENERGY USING THE PLASMON POLE 

APPROXIMATION 



After the angular integration which was performed in equation (|14|) . the imaginary part of the self-energy at the 
quasi-particle dispersion lu — e p reads 



Im£(p,p 2 ) 



4i f°° dq_ 

4f> Jo 



Erf 



Erf 



q 2 + 2pq + w q 

V^CTp 

q 2 + 2pq - L0 q 



-Erf 



Erf 



q 2 - 2pq + ujg 
q 2 - 2pq - o; q 

V2<T P 



n B (uJq) exp(w q /T) 

n B (-w q ) exp(- Wq /T) I . (Al) 



This equation represents a self-consistent equation for ImE(p,w — p 2 ) — yj2/ na p . 

Our aim is to derive an analytic expression, that approximates the numerical solution of equation (|A1[) for arbitrary 
p. To this end, we first look at the case of large momenta, p>K, and later combine that result with known expressions 
for the limit of vanishing momentum p — > 0, to produce an interpolation ("Pade") formula that covers the complete 
p-range. 

We perform a sequence of approximations to the integral in (|A1|) . First, we observe, that at large p, the term 2pq 
dominates in the argument of the error function. We rewrite equation (|A1|1 as 



dq 



Erf 



o ?w q 

' 2pq \ 
V2a p ) 



7T W pl 



dq 

2 4p J quj q 



Erf 



- Erf 



2Erf 



2pq_ 
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' -2pq \ 
V2a p J 

2pq ' 



Erf 



-2pq 



"bK) exp(w q /T) 



\/2<7 v 



K V2<J P 
nB(-w q ) exp(-cj q /T) | 

KW exp(w q /T) - n B (-u; q ) exp(-w q /T)] 



(A2) 
(A3) 
(A4) 



The integrand in equation (|A4j) contains a steeply rising part at q < —a p jp and a smoothly decaying part for at 
large q, i.e. when q ^> —<r p /p. Therefore, we separate the integral in equation into two parts, one going from q = 
to q = q = —<J p /p and the other from q to infinity. In the first part of the integral, the values for q are so small, that 
we can replace the plasmon dispersion by the plasma frequency uj p \. In the second term, the argument of the error 
function is large and the error function can be replaced by its asymptotic value at infinity, lim^^oo Erf(x) = 1. This 
leads to 



dq 



2 4p \ J quj pi 



2Erf 



_2pq L 
V2~p p 



[n B (w p i) exp(w p i/T) - n B (-w p i) exp(-w p i/T)] 
dq 



quj q 



KK) exp(w q /T) - n B (-w q ) exp(-w q /T)] 



(A5) 



Finally, we expand last term in powers of w q /T, which is justified at low densities (w q tx w p i), and keep only the first 
order, 



2T 

«B(w q ) exp(w q /T) - n B (-w q ) exp(-w q /T) = h 0(wq) 



We obtain 



dq 



2 4p 1 J qu! p \ 



2Erf 



2pq_ 
\p2o p 



[n B (w p i) exp(w p i/T) - n B (-Wpi) exp(-£d p i/T)] + AT 



dq 

9 W q 



(A6) 



(A7) 



Both integrals can be performed analytically, 

r*Erf =-2^ 2 F 2 (l/2,l/2;3/2,3/2;-2pV/^) 



2 . (A8) 



2d - 2 F 2 (l/2, 1/2; 3/2, 3/2; -2) = -1.3357 , 



where g = —a p jp was used. Note that in the second integral, the g term in the plasmon dispersion Q is omitted. 
2^2(01, 0-2', b\, 62; z) is the generalized hypergeometric function [461 ]. 
We arrive at the equation 

a p = _1 - 3357 Y^"^ I n B(w P i) exp(w pl /T) - n B (-w p i) exp(-w p i/T)] - y^^M 1 + K V/^) • (A9) 

At large p, the term K 2 p 2 /a p dominates the argument of the logarithm, i.e. we can write ln(l+« 2 p 2 /a 2 ) ~ ln(re 2 p 2 /cr 2 ). 
Then, we arrive at equation (|Alj) . given in section llVl 



APPENDIX B: PADE APPROXIMATION 



From the knowledge of the behaviour of a p in the limits p — > and p — > 00, a two point Pade interpolation formula 
can be constructed. For the value of the quasi-particle damping width at p — we take the expression 



on = v kT , 

2 

which is the exact solution of the self-consistent Born approximation (35| . 

The Pade interpolation formula is constructed in the following way: We make the ansatz 



(Bl) 
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p 1 + hp + b 2 p 2 

where the function <p(p) contains the logarithmic terms present in the behaviour of a p at large p, c.f. equation ([T 



(B2) 



~ , - ln£ ln£ ln£ 31n 2 £ ln 2 £ 

£ - ln£ + - -A + + — =A 

£ £ 2 £ 2£ 3 2£ 2 



In 3 1 



£ = ln(e + d -up 2 exp(A/T)/T) . 
The coefficients ap, ai, 61, 62 arc determined by power expansion at p — and p — » 00, 
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(B3) 
(B4) 

(B5) 
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and comparison to the behaviour of a p in these limiting cases, e.g. equation (|16p for large p and equation (jBljl for 
p — > 0. Setting the slope of er p at p = to zero, as well as the coefficient in front of the p~ 2 -term of the asymptotic 
expansion, we arrive at the following equations for the coefficients of the interpolation formula, 



a = -— vkT 



Ol 



a\ — aob q = 
a b 2 — a\bi = 



(B7) 
(B8) 



The solution reads 




a± = — k 




(B9) 




(BIO) 



which is given as equation Q17p in the main text, section (|TT| . 
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